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ABSTRACT. This paper details a description of the pattern of galaxy image distortion over the entire sky caused by the 
gravitational lensing which is the result of large scale inhomogeneities in our universe. We present a tensor spherical 
harmonic formalism to describe this pattern, giving many useful formulae. This is applied to density inhomogeneities, 
where we compute the angular power spectrum of the shear pattern, as well as the noise properties due to finite galaxy 
sampling and cosmic variance. We show that a detectable level of shear is present for very nearby galaxies, z < 0.2. 

, For such a shallow sample much of the largest signal-to-noise comes from very large angular scales, 6 > 10°, although 

' it is in the form of very small shear at a level < 10~ 3 . 
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1. Introduction 



In recent years gravitational lenses have provided a set of extremely useful tools for understanding the universe 
_l , around us. One of these tools, sometimes called the orientation correlation function or OCF, is a method whereby one 
£> ' searches for alignments in the orientations of galaxies on the sky (Tyson, Valdes, and Wenk 1990). Such alignments 
ON , will be caused by the deflection of light by the gravitational field of mass concentrations in front of the galaxies one is 
observing. In the weak lensing approximation the OCF is determined by the shear of the image deformation caused 
by this bending of light. Large format CCD's have made such observations of weak lensing shear possible and there 
' have been numerous studies of the mass concentration in clusters of galaxies with this technique (Tyson, Valdes, and 
Wenk 1990, Bonnett, Meiler & Fort 1994, Fahlman et al. 1994, Smail, Ellis & Fitchett 1994, Smail & Dickinson 1995, 
ON ' Tyson and Fisher 1995, Squires et al. 1996a&b, Luppino & Kaiser 1996). It is in the direction of clusters of galaxies 
at moderate redshift where the shear is liable to be greatest and the size of the central parts of such clusters is a 
good fit with the typical field of view of many CCD cameras, i.e. about 5'. Away from the direction of clusters of 
q , galaxies the shear is liable to be much smaller. In principle one should be able to detect very small shears by looking 
' for alignments among a greater number of background galaxies. One might obtain these greater numbers either by 
\ taking deeper images, which has the advantage that shear will increase with depth, or by looking at larger areas on 
■ the sky. With the advent of very large CCD mosaic cameras with a large field of view one can expect much wider 
, area surveys looking for weak lensing shear. A limiting factor in such wide surveys is the small amplitude of the shear 
one is trying to measure, and further study is required to see just how small one can reduce systematic errors when 
looking for galaxy alignments. 



One of the most interesting applications of weak lensing which is to map the mass distribution as close to the 
local neighborhood as possible. This will allow a comparison of the galaxy distribution and the mass distribution 
which should help us to understand galaxy formation and biasing. Gould and Villumsen (1994) have pointed out that 
the Sloan Digital Sky Survey (see Kent 1994 for a description of the "SDSS") which is imaging one quarter of the sky 
should be able to measure the mass distribution around the Coma cluster. Whether this goal is achievable depends 
on the level of non-correctable systematic errors, however even if present imaging surveys are not successful one can 
expect that very large area weak lensing surveys will be achieved sometime in the future. 

Predictions of the angular distribution of shear and amplification on the sky has been made in the context of 
the small areas on the sky where the small angle approximation is utilized (Blandford et al. 1991, Miralda-Escude 
1991, Kaiser 1992, Kaiser and Squires 1993). This is quite natural given that most weak lensing observations would 
be limited to extremely small patches on the sky. However looking ahead to a time when the area on the sky surveyed 
becomes large, such as in the SDSS, one will need to go beyond the small angle approximation in order to characterize 
the weak lensing pattern one observes. This is not to say that the lensing deflection angles will be large, they will not, 
but rather that when the area on the sky becomes large one needs to take into account the "curvature" of the celestial 
sphere when computing things like shear-shear correlation functions. Since the shear is not a scalar quantity, but 
rather a rank-2 tensor, one can not simply use the the scalar spherical harmonic expansion. In this paper is presented 
a tensor spherical harmonic expansion which the author feels is rather well suited to describing the shear pattern. 
Most published expositions of tensor spherical harmonics were motivated by describing 3-d tensor gravitational fields 
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and these expositions are complicated by the additional formalism needed to describe the 3rd dimension (see Thornc 
(1980) for a review). The formalism described here is rather simpler as it is restricted to symmetric traceless tensor 
fields on the 2-d sphere. Although developed independently this formalism is very close to that of Zerilli (1970). 

The paper is arranged as follows. In §2 we describe the decomposition of the shear pattern into two geometrically 
distinct types in a formalism which is easily generalizable from the small-angle approximation to the celestial sphere. 
In §3 is considered the various approximation used when computing the image distortion in lensing. In §4 is presented 
the tensor spherical harmonic decomposition of the shear, with a large number of useful identities. In §5 the tensor 
harmonic expansion is used to work out the formalism to describe the shear pattern produced by density perturbations 
in an Einstein-deSitter universe. In §7 we estimate the accuracy one might expect to obtain in measurements of the 
shear due to the "shot noise" from the finite number of galaxies that are available to probe the shear, and to "cosmic 
variance" (finite sampling). In §7 the formalism developed is applied to a phenomenological model of the density 
perturbations in our universe and illustrate the type of shear pattern one might expect to find and how it compares 
to the shot noise and cosmic variance. In §8 we summarize the results. 

2. Scalar and Pseudo-Scalar Shear 

The shear gives the degree to which images of objects we see on the sky are elongated, after factoring out the 
intrinsic shape of the actual objects. The shear at a point has an amplitude and a direction and may be described by 
a symmetric traceless matrix, i.e. 

( cos 2^ sin 2^ \ 

lab = 7 • o o • t 2 - 1 ) 

\ sin 2tp — cos 2<p J 

Here 7 measures the amplitude and tp the direction wrt to some fiducial position angle on the sky. Of course the 
direction is modulo n not 2ir since stretching an image in one direction is the same as stretching it in the opposite 
direction. Note that rotating the direction of the shear by 90° changes the sign of the the shear, while rotating the 
shear by ±45° produces a shear which is orthogonal to the original, in the sense that r y ab, y' ab = 0. (The Einstein 
summation convention is used here and throughout this paper.) Whether one rotates 45° to the right or 45° to the 
left changes the sign of the shear. One can construct any shear matrix by linear combinations of a given shear matrix 
and one rotated by 45°. 

Here we are concerned with shear fields on the sky. In this section we consider the small angle approxima- 
tion where the sky is approximated as a Euclidean plane. One may write any shear field on a plane as a Fourier 
decomposition, e.g. 

MZ) = e** (2.2) 
where q is the 2-d wavenumbcr on the (planar) sky and 

fab = d^ £%= (-°i o)- (2 ' 3) 

The tensor e a b is known as the Levi-Civita symbol in 2-d. The 7® modes have shear directed parallel to q while the 
7® modes have shear rotated 45° to the right from the q direction. This is all one needs to construct an arbitrary 
shear field. 

Since the 7® modes are just given by derivatives of the scalar mode function f(q,x) we will refer to this part 
of the shear as the scalar shear. The 7® modes are a geometrically distinct component of the shear, given by the 2nd 
derivative matrix of / but multiplied by the e c b which has the effect of rotating the scalar shear by 45° to the right in 
a right-handed coordinate system or 45° to the left in a left-handed coordinate system. A change in the handedness 
would have the effect of multiplying the pseudo-scalar shear tensor by —1. This sign difference depending on the 
handedness of the coordinate system means that the component of the shear transforms as a pseudo-scalar rather 
than a scalar, and it is called the pseudo-scalar component of the shear. The two types of shear pattern are illustrated 
in fig 1. One can represent any pattern of shear by a sum of scalar and pseudo-scalar shear 

The decomposition we have just described is similar to the decomposition of a vector field into it's vortical 
and non-vortical parts. Note however that in 3-d for each Fourier mode their are two linearly independent vortical 
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components and one must choose a azimuthal angle about the wavenumber to specify one. For this reason one cannot 
write these vortical modes in terms of derivatives of the scalar modes e lk x plus a handedness as this will not specify 
the azimuthal angle. In 2-d one may express arbitrary vector and tensor fields in terms of scalar functions Levi-Civita 
symbols, as we have done above. There is not the distinction between scalar, vector, and tensor components of fields 
higher dimensions. 

In the usual treatment of gravitational (or non-gravitational) lcnsing, one thinks of the lcnsing as being a 
displacement field on the sky, i.e. in a given direction on the sky what one sees is at a true position in space which is 
the apparent position displaced by some amount, A. This is what we will later call the displacement approximation. 
For arbitrary gravitational perturbations it is problematic to define what the apparent position is a hence what the 
displacement field is. In the next section it is argued that for weak fields it is an excellent approximation to treat the 
lcnsing as a displacement field on an unperturbed space-time and this is the approximation we will use here. Let us 
denote the displacement from apparent to true position on the planar sky by A. The 2x2 deformation matrix is 
given by 

tpab = Sab + A a , fc = (1 - k) 8 ab - 7ah + u e ab (2.4) 

which we have decomposed into term which gives the trace, a term which gives the symmetric traceless part, and a 
term which gives the anti-symmetric part. The trace gives the expansion, k, the traceless symmetric part give the 
shear, j ab , and the anti-symmetric part gives the rotation, ui. The image amplification is 1/((1 — k) 2 — 7 2 + w 2 ) or 
approximately 1 + 2k in the weak lensing limit. An arbitrary displacement field can also be written as a Fourier 
integral of scalar and pseudo-scalar components: 

A B (f) = ~ J d 2 q [<f>®(q) f, a + cf>®(q) f, c e c a ] f(q, x) = e^ S . (2.5) 

Note that f. c e c a rotates the gradient of / by 90°. When the deformation is given by a displacement field on a Euclidean 
plane we find the following relations 

V 2 k = +1 a \ ab V 2 uj = -7 afc ,a C e c fc. (2.6) 

Since V 2 is invertible the scalar shear is given by the expansion, n, while the pseudo-scalar shear by the rotation, oj. 
A scalar and pseudo-scalar displacement field is illustrated in figure 1. 

It is unlikely that anything other than density inhomogeneities contribute significantly to the shear in our 
universe. These density inhomogeneities produce scalar metric perturbations. One shouldn't confuse the "scalar" -ness 
of density perturbations which has to do with transformation properties in 3 spatial dimensions, with the "scalar" -ness 
of the shear which has to do with transformation properties on the 2-dimensional sky. Nevertheless it is true that 3-d 
scalar perturbations will only produce 2-d scalar shear, at least in the weak lensing limit. In contrast tensor and vector 
metric perturbations will contribute to both scalar and pseudo-scalar shear, however we expect this to be very small. 
Thus we really do not expect there to be any significant pseudo-scalar shear. We have chosen to include pseudo-scalar 
shear in our analysis for a variety of reasons. Most importantly because the pseudo-scalar shear is part of what one 
measures when one looks at galaxy alignments, and one should not leave it out of ones analysis. Also, one never 
knows, maybe the pseudo-scalar shear field is not negligible, one should measure it and see! 

Kaiser (1992) noted that the redundant information in the shear field could provide a useful check of ones 
observation. One can determine the pseudo-scalar u> in the same way one determines the scalar k if one first rotates 
the galaxy position angles by 45°, and then one should check that the derived w is is consistent with zero. (Kaiser et 
al. 1994; A. Tyson private communications; Stebbins, McKay, and Frieman 1996). Just such a procedure has been 
implemented in Luppino and Kaiser (1996) where they compare j ab , a b to j ab ^ ac e c b showing that the latter is indeed 
much smaller than the former in the field of the cluster of galaxies msl054-03, and consistent with zero. 

3. Approximations 

In this paper we make predictions of extremely small image deformation. These are calculated using a number 
of approximations and one must be careful that the approximations are accurate enough so that these predictions are 
accurate. In other words one must be sure that there are not other small but significant contribution to the lensing 
which are not included. Many approximations become more accurate as the image deformation becomes smaller, while 
others do not. In this section wc identify most of the approximations used here. None of these approximations are new 
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to this work, and not all studies of gravitational lensing make use of all of the approximations used here. Below we 
will use the geometric optics approximation in weak field gravity for perturbations from Friedman-Robertson- Walker 
space-times. These are clearly excellent approximations for cosmology. It is generally assumed that we are in the "weak 
lensing" limit, namely that the deformation tensor of eqs (2.4) & (4.1) are close to the identity matrix, although many 
of the formulae are true even for strong lensing. For isolated regions around galaxy clusters, galaxies, and stars weak 
lensing will not be a valid approximation, and generally as one goes to smaller angular scales the stronger the lensing 
becomes. To what extent effects from small scale strong lensing creep into larger scales is a subject which needs further 
study. 

When computing image distortion we integrate the perturbation to a beam of photons along the unperturbed 
path of the photons. To do otherwise would be to include terms which are formally 2nd order in the metric perturbation 
so one can argue that this is part of the weak lensing approximation. Even when the image distortion is large this may 
also be a good approximation if the statistical properties of matter field which one traverses on the perturbed trajectory 
are quite similar to that traversed on the unperturbed trajectory. However the argument that density perturbations 
only produce scalar patterns of shear rests largely on the assumption that it is valid to use the unperturbed trajectory. 
In general for strong lensing one can expect density perturbations to produce both scalar and pseudo-scalar shear. 
An exception occurs when the mass which produce the strong lensing is localized to a narrow range of distances. In 
this case one can still use the unperturbed trajectory to compute the total deflection accurately. Thus in studies 
of lensing around individual galaxies or galaxy clusters one can assume that the shear pattern is purely scalar since 
it is unlikely that two large mass concentrations will line-up sufficiently to contributed significantly to the lensing. 
However on small enough angular scales the strong lensing will be due to a variety of objects along the line-of-sight 
and the shear pattern will contain both scalar and pseudo-scalar components. Again further study should be done on 
the contamination by small scale strong lensing of larger scale weak lensing. 

One approximation which is important for this paper is what will be called the deflection approximation. 
Namely that one can think of the effect of gravitational lensing as a deflection, where the light has been deflected 
from a "straight path" in an underlying flat geometry. While this is a perfectly valid description for lenses made of 
refractive materials in flat space, it is only an approximation for gravitational lensing. This is because there really 
is no underlying flat geometry, after all the light is deflected because the space is not flat. Below we will argue that 
in the cases of interest the space is flat enough so that this deflection approximation is a good one. Even though in 
this paper we consider very weak lensing, the image distortions we consider are still much weaker than the intrinsic 
distortion pattern caused by the curvature of spacctimc. Note that this is not an intrinsic property of the weak lensing 
approximation, but rather has more to do with how the amplitude of metric perturbations in our universe vary with 
scale. If we had very large perturbations entering the horizon today much of the image distortion would be a result of 
the global geometry rather than the deflection. 

Let us try to make these ideas more precise. In §5 of this paper we consider the lensing due to density 
perturbations. The deflection as a function of distance from the observer is computed using the metric evaluated 
on the unperturbed path, i.e. by twice integrating the "gravitational acceleration" of the photon perpendicular to 
the unperturbed path. With this deflection we compute the intersection of both the unperturbed trajectory and the 
"deflected trajectory" with a surface defined by a fixed coordinate distance from the observer, i.e. r = \J x 2 + y 2 + z 2 , 
which we will call the "source sphere" . We use the the Jacobian of the mapping of the Newtonian coordinates from 
the deflected position on the source sphere to the undeflected position as the deformation tensor in eqs (2.4) or (4.1). 
This procedure has some obvious failings. Firstly if the photons do not intersect the source sphere at right angles then 
even for photons which are not focused or sheared the Jacobian matrix will contain an eigenvalue of 1/ cos 9 where 8 
is angle between the normal to the surface and the direction of the perturbed photon trajectory. Thus even if there 
is no true image distortion one would infer a de-amplification and shear. This error is not large, in the sense that 
the leading error is cx 8 2 , so is an insignificant error in the sense of weak lensing. In cosmology the deflections are 
rarely greater than 1CP 4 radians even when the lensing is strong. A more worrisome error comes from the fact that 
the source sphere defined above is not truly spherical in geometry. Scalar perturbations in Newtonian gauge have an 
isotropic metric so the intrinsic geometry of the source sphere does not lead to any spurious shear. However a small 
coordinate patch spanning a coordinate area A on the source will actually spans a physical area (1 — 2$/c 2 )A where 
$ is the Newtonian gravitational potential. Unless one corrects for this one will add a spurious contribution to the 
expansion, n, of 2$/c 2 . Since this contribution to k is linear in the metric perturbation it does not become arbitrarily 
smaller than the true signal in the weak lensing limit. However, for our own universe we know that $ < 10~ 5 c 2 while, 
as we will see below, k > 10 . Thus by ignoring the intrinsic geometry of the source sphere we only make a small 
error. As mentioned above, the reason for this fortuitous numerology has to do with the shape of the power spectrum. 
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No such fortuitous numerology applies to the case of gravitational lcnsing from gravitational radiation, so one cannot 
use the deflection approximation in this case. We do not consider lensing from gravitational radiation in this paper. 
We will ignore the intrinsic geometry of the source sphere in the rest of the paper. While it is not terribly cumbersome 
to include this term, it would invalidate the simple and accurate relationship between the different components of the 
image distortion, k, to, and "f a b as expressed in eq (2.6), or implicitly in eq (4.10). 

A more general treatment of image distortion, can be gotten from the optical scalar equation (e.g. Hawking & 
Ellis 1973; Schneider, Ehlers, & Falco 1992) whereby one measure the distortion of a bundle of light rays with respect 
to an orthonormal basis parallel transported along with the light rays from the observer. This approach makes no 
approximation other than assuming the validity of geometric optics. It would be too cumbersome for our purposes 
to use this approach on unperturbed trajectories, i.e. to use ray-tracing. One could approximate the optical scalar 
equation by implementing the same procedure along an unperturbed trajectory. One property of the optical scalar 
equation which is something of a drawback is that since the parallel transported basis vectors will rotate with the 
photons, this method does not yield any information on image rotation. While in the case of strong gravitational fields 
it is difficult to make a sensible definition of image rotation, in the case of weak fields where the image distortion is 
generally much larger than the metric perturbation there is a real sense in which images may be rotated and it would 
be useful to include this effect in ones computations. 

In §5 and afterwards a set of assumptions and simplifications are implemented. We do only consider a flat 
matter-dominated Einstein-deSitter cosmology, although one could clearly generalize all of formulae to open or closed 
cosmologies with or without a cosmological constant. The formulae used in §5&§7 assume that the perturbations grow 
according to linear theory. This is a good approximation on large scales > 10/i _1 Mpc but we when we consider the 
shear at shallow depths we do venture into a regime where linear theory is not appropriate. In fact the model power 
spectrum in §7 is no-linear on scales of interest. One could try to correct this by having the power spectrum vary as 
one looks back in time, i.e. with radius. Since the results presented here are meant to be illustrative this approach 
is not warranted. In any case since the emphasis is of this paper is on shallow surveys these sorts of corrections are 
liable to be small. 

4. Decomposition of Shear on the Celestial Sphere 

Even though the coherence in the shear field falls off with separation there will be some small coherence even 
over very large angles on the sky. To understand this large angle coherence one must take into account that the sky 
is a sphere and not a plane. Here we consider the weak lcnsing on the celestial sphere, i.e. the sphere of the sky. 
For the moment we consider the image deformation at a fixed distance but later generalize the analysis to include the 
fact that the galaxies whose images are distorted span a range of distances. Thus we have the a mapping from the 
direction in which you are looking, n, to the direction you are looking at, m, at a given distance, r. The effect of the 
lcnsing is to displace the photon trajectory at a given distance by the 3-d vector A = m — n. If the angle by which 
the light rays are bent are sufficiently small then A and n are nearly perpendicular so that we may approximate A 
as being in the tangent space of the the direction sphere h, i.e. A(r) is a vector field on the sphere. As in §2 we may 
define the deformation tensor and it's various components 

g a b + A a:fc = ipab = (1 - k) g a b - "fab + ue a b (4.1) 

where and A a: b refers to a covariant derivative of A wrt the metric g a b on the sphere. 1 Here k, j a b, and w give the 
expansion, shear, and rotation. The shear is defined such that g ab "f a b = which fully specifies k. From the shear field 
one may define two new quantities 

7s = V- 2 7ah :afc 7P - V- 2 7Q6 :fcc e a c . (4.2) 
where V~ 2 is the inverse Laplace operator on the sphere, i.e. 

V 2 cj> = cj>;a a V~ 2 (f>;a a = 4> + constant. (4.3) 

We must add the arbitrary constant because the V -2 is not completely defined since V 2 has a zero eigenvalue which 
corresponds to the constant, I = 0, eigenmode. Here I gives the eigenvalue of the Laplace operator which is —1(1 + 1). 

1 All of the indexed quantities are components of tensors on the sphere. Below we will use the coordinate bases defined by spherical 
polar coordinates, (8, <f>) on the sphere. Although the two basis vector are orthogonal they are not orthonormal, and in particular, tensor 
components with a <fi index will contain extra factors of sin 2 6 when compared to tensors defined with respect to an orthonormal basis in 
the same direction. 



5 



The quantity 7 S is contributed to only by the scalar part of the shear, as described in §2, while the quantity 7 p is 
contributed only to by the pseudo-scalar part of the shear. From 7 S and 7 P one can fully reconstruct jab- However 
from 7 a h or from 7 S and 7 P , one cannot, in general, fully reconstruct the displacement vector A. As we shall see, the 
I = 1 eigenmodes do effect A but not -f ao . The fact that the I = cigenmodes also do not contribute to j ao does not 
matter since it also does not contribute to A. By definition n and u have zero mean and therefore no contribution 
from the I = eigenmode. It therefore makes since to define 7 S and 7 P to have zero mean and therefore no contribution 
for the I = eigenmode. By construction they have not contribution from the I = 1 eigenmodes. 

Above we have used the completely antisymmetric tensor defined by 

which commutes with covariant differentiation, i.e. e ab :c = 0. When e ao is contracted with a 1-indcx vector it has 
the effect of rotating the vector by 90°. Whether the rotation is to the left of right can be chosen arbitrarily, but 
whatever choice one makes the sign of e ao is multiplied by —1 under a parity transformation, and thus a leftward 
rotation is transformed to a rightward and vice- versa. Thus quantities containing odd numbers of e b's are odd under 
parity and we refer to them a pseudo-scalar part of the shear. Note also that the shear e a 2 b^ a b yields another traceless 
symmetric tensor which if interpreted as a shear tensor has it's shear rotated by 45° wrt that of the original shear 
tensors. Rotating the displacement vectors by 90° corresponds to rotating the shear by 45°. Performing either of these 
rotations will transform a purely scalar patter to a purely pseudo-scalar one, and vice-versa. 

Spherical Harmonic Expansion 

In analogy with the Fourier expansion of §2 we may expand the displacement vector on the sphere in terms of 
a spherical harmonic functions Y (!im) (n). These functions are defined such that 

V 2 y (; , m) (n) = -i(/ + l)y (I , m) (n) j d 2 hY^ m) Y;, tml) = 6 w 5 mm , 

Using 2-d spherical geometry this implies the identities 

Y {l , m ;- ab .. ab = v 2 v 2 y (; , m) + v 2 r (1 , m) = i(i + i) (i(i + 1) - i) y { 

J d 2 hY {l , m):ab Y i ;, m ,- ab = 1(1 + 1) (1(1 + 1) - 1) 5 w 6 mm , 
The spherical harmonic analog of eq (2.5) is 

oo + 1 

A « = - E E «™> y «->- + < m) y^y. b e h a ) (4.7) 

l—l m— — l 

including both scalar and pseudo-scalar displacement. Note that we do not include a I = m = term in the sum since 
Y(o,o):a = 0. Furthermore since A a is real we require 

<t>% m) = (-ir< m) * $?. m) = (-ir<^L>* smce y*_ ra) = (-iry (I . m) . (4.8) 

In the above decomposition the 0® m) terms give the scalar part of the displacement vector which is proportional to 
the gradient of the spherical harmonic, Y (hm) . The pseudo-scalar part, given by the <j)® terms, is proportional to the 
gradient of F (!im) rotated by 90°. The deformation tensor is then given by 

oo +Z 

i>ab = gab (<l>l m) Y V,my.ab + <, m) ^ (! , m y.ac^ ''b) (4-9) 

/ — 1 m— — l 
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(4.5) 



(4.6) 



and it's component parts by 



«=4ee i{i+i)4>t m) Y {l , m) 

l — l m— — l 

w = - 2^212 Kl + ^l m) Y {l , m) (4.10) 

l — l m=—l 

°° +l ( 1 1 \ 

lab = E E ( ^®,™)( y («.™)=«6- ^9abY {l:my . c - c ) + cj)f lm) -{Y {l:my . ac e% + Y {ltmy . bc e\)\ 



where we have used eq (4.5). In the shear we do not include the I = 1 terms, since the deformation tensor from these 
modes do not contribute to the shear as we will see below. Combining (4.2,6,&10) one finds 



00 +i 



7s = (-v- 2 </w a « + K ) = iE E a + 2)(/-i)< m) r (i 

1=2 m=-l 

7 P = -(V-Va h : " c 6 a c - V" 2 (a;e a6 ) :bc e a c ) =- £ E (* + 2)(i - 1) < m) ^> 



m) 



(4.11) 



= 1 m— — l 



and that the Z = 1 terms do not contribute to either 7 S or 7p . For large Z we recover the small angle result, i.e. 7s = k 
and 7 p = —ui. These relations hold for any deformation of our sky onto the 3-d space, but they are nevertheless 
they are non-trivial consistency checks. For example one can hope to measure both 7 S and k, the former from galaxy 
alignments and the later from galaxy brightness distributions in various bands. These being independent measures of 
the same thing would provide a consistency check. If they were found to be inconsistent then would discover that one 
is not measuring variations in alignments and brightnesses due to an image deformation. Measurement error, Galactic 
extinction, intrinsic alignments, intrinsic galaxy clusterings and other effects will all contribute to inconsistencies. 

Decomposition of Distortion Pattern 

Given combinations of A a , k, w, j s , 7p, or j ab one can compute the mode coefficients via 



1(1 



TT) / dHY ^- Aa = +/GTI) / d * AY ^ Aa - = -WTT) I d2hY ^ K 

^/^n^ ):o A^ = -^/^nr* m) A^ = ^/ 



d 2 nY ( l m) co 



= + /(/ + l)(/J + l)-l) S d ^n,^r b = J(n _\y- 1 J<PnY (lm) (7s + n) 



1(1 + 



I / d 2 hY* 

i)(i(i + i)-i)J dnY ^y- 



c ni.ab 



1(1 + 



^y-y / d'nY ( l m} ( 7p - w ) 



(4.12) 



+ 



(l + 2)(l + 



,) is 



(I + 2)(l + 1)1(1-1) 



cd _ ad 



d 2 nF ( * m):ac e^7 



(l + 2)(l 



— Jd 2 nY ( l m)lp 



using the orthonormality relation of eqs (4.5&6). We see that one cannot obtain the I = 1 terms from the shear since 
substituting I = 1 in the last four expressions one finds that the prefactors are infinite while the integrals are zero. 

Mean Square Expansion, Rotation, and Shear 

For a given realization of the shear pattern one can construct the quantities 

l 2 (l + lf — — 12(1^^2 • 



4(2/ + 1) 



m——l 



(4.13) 



m— — l 
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The notation is used to indicate that these quantities may be used as estimators of the power spectra coefficients, 
Cf and Cf, defined below. Combining the decomposition in eq (4.10) with the orthonormality relations of eqs (4.5&6) 
we find that the mean square k, oj, and j a b averaged over the sky is 

^ = E^^ (4-14) 

l=i 

where ||7 a h|| indicates the determinant of the matrix 7%. We use the notation — to indicate averages over the sky. One 
expects the mean square shear to receive it's largest contribution from small angular scale so that the large I terms 
dominate these sums. If this is the case then 7 2 w k 2 + u 2 . One could construct equations analogous to eq (4.14) for 
mean-square values after convolution with a window function by one could include a \ Wi\ 2 factor to each of the terms 
of the sum, where Wi is the spherical harmonic decomposition of the window function. 

Explicit Construction of the Shear Tensor 

Here we give the explicit expression for the shear tensor in terms of the spherical harmonic mode coefficients. 
Using spherical polar coordinates, (#,</>), the metric and anti-symmetric tensor are 

9ab= U sin 2 - 1 9 = S ' m2(> t% = I I !4 i ,; 




so (4.10) becomes 

00 +; 



lab = + E 

1=2 m=- 

oc -\-l 

+E E 



ll ' m> y Y (hm):Sl p \Y (l , m):< ^ - \ sin 2 6Y V , m) -.ee J 

& ( - ■±sY(i, m y.e<t> ±sin6%, m):ee - 3^ 

<i,m) Usin^y^,^-^^,^ sin 6 Y (ltmr . e<j> 



(4.16) 



Z=2m=-i \ 2 k..".;-™ 2sin 

To continue we need the explicit form of the spherical harmonics, which we take to be 



WM) = y^ i f T ^l P ™ (cos ' )eim0 (417) 

following the conventions of Jackson (1975). 

Using this explicit form of the spherical harmonics one can show that the deflection vectors Y {1}0 y.be b a , Q^(i,+i):6 + 
Y(i,-i):b)e b a , an d i(Y(i,+i)-.b— Y (1 _ 1): b)e b a , are generators of solid body rotations of the sphere. Clearly solid body rotation 
generates no amplification or shear, but does translate and/or rotate the images and thus we see why the I = 1 pseudo- 
scalar displacement contributes nothing to the shear, j a b- The I — 1 scalar displacement is just the pseudo-scalar 
displacement rotated by 90° which would yield a shear rotated by 45°. However since the I = 1 pseudo-scalar shear is 
zero, so must be the I = 1 scalar shear. 

Let us defining two auxiliary function via 



'2Z + l(Z-m)l , , „, lm<t> 1 1 v 

(4.18) 



l 2l + l(l-m)\ (cQsd)em ,_ Y 
4tt {l + m)r^ [C ° hd)e _y<i - ): 
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or more explicitly 



G^ m) {x) =m ((/ - 1) Pr{x) -(l + m) P^x) 

we find from substituting (4.18) into (4.16) that 



(4.19) 



in— — I \ 



760 _ 70fl 1 21 + I ^ /(I - m)\ 



(4.20) 



EVlT E J 77x3 (^^(cos^ + ^^G+^^cos^)) e 



sin6> sin# ^ V 47r ^ \/ (Z + m)! 

Z=2 m=-Z V V / 



Note the symmetry relations 



^-^(-irf^qL, (4-2i) 

and the orthonormality relations 



(4.22) 



From these identities one can rederive the formula for j 2 in cq (4.14). 

Below we will find that it is of particular interest to calculate the shear at the pole of the coordinate system. In 
spite of the fact that the coordinates system is singular at the pole, mode expansion at the pole takes a particularly 
simple form. The coordinate singularity means that one should consider only -fgg in the limit 9 — > for a given value 
of <fi. In this limit the value 700 will retain an 0-dependence and from this one has the angular dependence of the 
shear. One finds that in the 8 = limit that only m = ±2 contributes 



<# m >(l) = if^fj'Wi^,. , (4.23) 



which on can use to compute the 2-point correlation function of the shear. 

Correlation Functions and Power Spectrum 

Henceforth we shall explicitly include the radial r-dependence of the deformation since we may be interested in 
correlating the different components of the correlation at different redshifts. We shall also be considering a statistical 
distribution of deformation fields. Averaging over realizations, one can characterize any set of 2-point functions in 
terms of the power spectra which we choose to define by 

C>)^, m ,V)) = ^0S w 8 mm > <,„,>')) = ^0S w 8 mm >. (4.24) 

The Kronecker ^-functions and the lack of m-dependence follow from the assumed homogeneity and isotropy of the 
distributions we are considering. We include a dependence of the deformation with depth r. We will also assume that 
the distribution of deformations is even under parity as it would be for gravitationally induced deformations. Since 
4>® is even under parity and <j>® is odd under parity it follows that 



v < m >W< W) V))=0. (4.25) 
Both C®(r, r') and C®(r, r') are real and are non-negative when r = r' . 



The correlation functions of various quantities can be computed in terms of the correlation function 



1 oo 

C K (r, r',tf) = («(r, n) «(/, n')) = ^ £(2/ + 1) fl(cos0) Cf (r, r') 

oc 

a(r,r',i?) = (u;(r,n)u;(r',n')) = — £(2i + 1) P t (cos 0) Cf (r,/) 



47T ■ 

£=1 



(/c(r,n)w(r',n')) = 
2Z + 1 



(4.26) 



(C©(r,r')G+ 2) (cos^) + Cf (r,r')G- 2) (cosi?)) cos 2^ cos 2<p' 

+(Cf (r, r') G" 2) (cos 0) + Cf (r,r') G+ 2) (cos *#)) sin2<^ sin 2p' 
1 °° 27-1-1 

C 7K (r, = -_^_±_C'®(r,r')P I 2 (co S <?) cos2^ 

(=2 ^ ' 

1 00 27 -1-1 

C 7U (ry,tf, V ) = +-J2 W f T) Cf(r,r')P^cos^ sm2f 

In each of the above expressions $ = Z(n, n'), i.e. the angle between the two points. C 7 measures the correlation of 
"fab with itself, while C 1K and C 1UJ measures the correlation of "f a b with k and u>, respectively. The meaning of if and 
if' are described in fig 2. Setting Cf = in C 7 we obtain a correlation function similar to those in Kaiser (1993). 
The C ++ (a), C +X (a), and C x x (a) of that paper corresponds to setting ip = if' = a, if = if' + j = a, ip = if' = a — \ . 
Kaiser (1992) paper did not include pseudo-scalar shear since this is not contributed to by density perturbations. 

From the above correlation function one may calculate the the mean square quantities at a given distance 

<. 2 ) = C K (r, r,0) = £ ^ Cf (r, r) <c 2 ) = G.(r, r, 0) = £ ±±± Cf (r, r) 
(7 2 > = {\l ab la b ) = G 7 (r, r, rf, 0, 0) + G 7 (r, r, 4,1,1) = g + l (j + 2)0- 1) r) + cf (r , r)] 

Note that the factor (7 + 2)(7 — + 1)) = 1 + 0(7~ 2 ), and since most of the shear is expected to be generated at 

very small scales we find 

( 7 2 }«( K 2 ) + ( W 2 ). (4.28) 

Remember that although (4.14) is quite similar to (4.27) the former refers to averages over the sphere for a given 
realization, which could be a observed quantity, while the latter refers to an average over realization and is a purely 
theoretical construct. The reasons that the formulae are so similar is because, by construction, 

(Cf (r,r')> = Cf(r,r') (Cf (r,r')) = Cf{r,r'), (4.29) 

i.e Cf and Cf are unbiased estimators of Cf and Cf, respectively. 

Comparing (4.26) with the orthonormality relations, (4.5&6) as well as 

i J dxP^P^x) = 2^1*"' ( 43 °) 
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we see that if one knows the correlation functions one may calculate the power spectra using 



Cf (r, r') =2tt / dx Pi (x) C K (r, r' cos™ 1 x) 



(Z + 2)(Z - 1) J_ x 



dx I dip Pf (x) C JK (r, r' , cos 1 x, ip) cos 2(f) 



/l 2 /"l 
dx Pi(x) C u (r, r' cos™ 1 x) = — — — — — dx dip P?{x) G 7W (r, r', cos™ 1 x, </?) sin 2</> 
-l (' + 2)(t — 1) J-i Jo 

4 1 f 1 /" 27r /" 27r 

Gf(r,r') =- ^ + 2 - )2 (; _ x)2 J x dx J o di P J o VC 7 (r,r', cos™ 1 x,ip,ip')x 

(G+ 2) (x) cos2<p cos2ip' + G™ 2) (x) sin2< ) 5 sin2y>') 

dip / rf<^' C 7 (r, r' , cos™ 1 x, </?, y/) x 
Jo 

(G+ 2) (x) sin2(/? sin 2<p' + G™ 2) (x) cos2(^ cos2<^') 



C?(r,r>) 



4 1 

^ (Z + 2) 2 (Z-1) 2 



dx 



(4.31) 



Small-Angle Limit 



Now let us consider the small angle limit of the formulae just derived. Over small angles the surface of the 
sphere is approximately a plane, and the wavenumber, I, can be consider modulus of a 2-d vector wavenumber in this 
plane in a Fourier decomposition of the displacement if I 3> 1. This Fourier representation is context in which most 
analyses of image deformation has been analyzed, and the formulae derived above should approach this Fourier limit 
for I 3> 1 over regions much smaller than one radian. In particular the angular correlation function should approach 
the Fourier limit for I 3> 1 and (?Cl. In the small- angle large-! limit one can show 



G± m) (cos0) = (-l) m ^l e - m (J m - 2 (W) ± J m+2 (W)) + 0(j, j,9 2 ) 

pr (cose) = (-i)-|±^|j m ^)r- + o(j±- 2 , ^ 2 ) 

Thus in the small-angle approximation the correlation functions become 

C K (rS,0)Ks— dllJ {W)Cf{r,r') G 7K (r, r', 0, ip) « -— / dllJ 2 (M) Cf{r,r') cos2ip 
C u (ry,4)si—J dllJ {W)Cf{r,r') G 7W (r, r', 0, w +— / tH / J 2 (W) Cf(r. r') sin 2^ 



(4.32) 



1 /" c 

c 7 (ry,i?, v ,/)«— y 



dZZ 



(4.33) 



(G®(r,r') (J (W) + J 4 (h?)) + Cf (r,r') (J («) - J 4 (W))) cos 2^ cos 2^ 

+ (G®(r,r') (J (W) - J 4 («)) + Cf(r,r') (J («) + J 4 (/<?))) sin2<^ sin 2^' 

while the small angle limit of (4.31) is 

Cf{r,r') w27T / J (lti)C K (r,r'd) w -2 / di?W dip J 2 {W)C 1K {r,r' ,$,ip) cos 2(j> 

Jo Jo Jo 

f-CO f-CO t>2-K 

Cf(r,r')^2n d-d J o (Z0) C u (r, r'ti) w -2 / d00/ J 2 (W) G 7W (r, r', 0, 93) sin 20 
Jo Jo Jo 

1 />oo />27r />27r 

G, ffi (r,r')~- / dtft? / dp/ G 7 (r, r', d,ip,ip') (J (W) cos 2(p - p') + ^(^) cos 2( ¥ > + y/)) 

~ Jo Jo Jo 

Gf(r,r')«- / di?7? / <W G 7 (r, r', 0, <p, (J (lti) cos 2(cp - 99') - J 4 (Z0) cos 2(y> + 95')) 

~ Jo Jo Jo 

Using a 2-d definition of a power spectrum where 

/•OO 

</(n) /(n')> = 2tt / dllP f (l)Jo(W) 0=Z(n,n') 
Jo 



(4.34) 



(4.35) 
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we see that 

P «(0 = 4^Cf ^(0 = ^Cf (4.36) 
and that 7 a fc is described by these same two power spectra, e.g. 

dll(P K (l)+P u (l))^—J dll{Cf{l) + Cf). (4.37) 

Visibility Functions 

To measure the shear one typically looks for correlations in galaxy ellipticity orientation, and one needs enough 
galaxies to find the correlation in the noise generated by the random orientation of the galaxies. Usually one doesn't 
know the distance to the galaxies except in a statistical sense, and even if one did, one would have to sum over galaxies 
at different distances to obtain a significant signal. Thus the shear at a given distance is not what one measures, but 
rather some weighted average of the shear at different distances. We may represent this distribution of distances by a 
visibility function, V(r), normalized so that 

pOO 

drV(r) = l, (4.38) 



/ 

Jo 



10 

so that the average expansion, rotation, and shear is 



P OO pOO P OO 

7t(n) = / drV(r) K (r,h) uJ(n) = / dr V(r) u(r, n) j ab (h) = / dr V(r) j ab (r, n). (4.39) 
Jo Jo Jo 

One can apply all of the formulae derived above by replacing the quantities at fixed distance with the same quantities, 
averaged over distance, e.g. the mode coefficients 



(!,m) 



roc _ poo 

= drV{r)4>l m) {r) 0® m) = / dr V (r) 0® m) (r) , (4.40) 
Jo Jo 



or the correlation functions 



where x might be k, w, or 7. 



pOO pOO 

CW(0, ■■■)=/ drV(r) dr' V(r') C x {r,r' ,0, • • •)■ (4.41) 
Jo Jo 



It is not the purpose of this paper to be terribly realistic about how one samples the shear with a given galaxy 
sample. The only visibility function which we will consider explicitly in this paper corresponds to uniform sampling 
in space up to a given distance, r max . The visibility function in this case is 

V(r) = 3^— (4.42) 

where r max gives the maximum comoving distance which one includes in ones sample. 

5. Application To Density Inhomogeneities 

Here we consider the image distortion from a linear growing mode density field in a matter-dominated Einstein- 
deSitter cosmology. For these growing modes the Newtonian gravitational potential is constant with time at a given 
comoving position. It is convenient to use the Spherical Harmonic - Spherical Bessel function expansion so the 
gravitational potential may be expanded 



00 1 



1 — ' L POO 

1=0 m =-l ^° 



*«, m ) (k) = ~ dr r 2 j t (kr) / d6 sin / d<f> Y* (0, 0) $(r, 0, 

1* Jo Jo J -it 



(5.1) 



0) 
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where r is the comoving distance from the observer. For the spatially flat geometry we are considering, these modes 
are orthogonal cigenfunctions of the Laplace operator. For an statistically homogeneous and isotropic distribution the 
correlation function of mode amplitudes is of the form 

(* ( ,. m) (fc) 5* 1 , im , ) (fc / )) = (47r) 2 P$(fc) 4 S(k - k') S w 5 mm , (5.2) 



where P<j> is the power spectrum of the potential fluctuations. Using the relation between the potential and overdensity, 

3 If? 
2 a 2 

we may also relate the potential power spectrum to the density power spectrum: 



V 2 $ = A-Ka 2 G5p = ^5, (5.3) 



9 (H^ 4 



Below we will use the notation n for (0,cp). Clearly these scalar modes will lead only to scalar components to the 
image deformation. Thus we may describe the deformation by a scalar function 

oo 

Aa = -0 :a = £ C„> W*) $L> = °" (5.5) 

1=1 

Comparing eq (5.5) with eq (4.12) we find 

e (r,n)-2 T° bS drl T - {v ° hs - p $( Xobs + ( Vohs - rj) n, rf) = 2 f dr' T —^- 3>(r', n), (5.6) 

Jrjobs-r r Wobs ?7 J JO rr 

and comparing this with (5.1) we see that 

f°° ~ r r - r' f°° ~ 

< m) (r) = 2/ dkk 2 $ (l , m) (k) dr' -j l (kr')=2 dkkH (l ^(k) I t (kr) (5.7) 

Jo Jo rr Jo 

where 

4(*) = / -(l-y)ii(^) «>o. (5.8) 

There is no need for an I = term since there is no Z = contribution to shear or expansion. In the appendix wc 
present analytic expressions for Ii(x) as well as series and asymptotic expansions. We will argue below that in most 
practical applications is well approximated by eq (A6), i.e. 

!«(*)« (at--) *(*--) o, = — ^ & ' = — T(|y- ( 5 - 9 ) 

which we call the asymptotic approximation. Combining (4.24), (5.2), (5.4), and (5.7) we find 

Cf(r,r') = 16ir 2 l 2 (l + l) 2 / dkk 2 J, (At) /;(At') P«(fc) 

Jo 

= 36tt 2 / 2 (Z + l) 2 g J,(fcr) W) P tf (fc) 



(5.10) 



HoY f°° dk ,- 



Taking into account a visibility function one finds 

Cf = 36n 2 l 2 (l + lf\^) j( ^[T^l 2 P s (k) (5.11) 
where 

f-OO 

h(k) = / drV{r)h{kr). (5.12) 
Jo 

From the asymptotic series of eqs (A4&5) we see that 



^ f 'W = (2TTW+i) ri "io drV{r)r & Il( * ) = & Il(fc) = Vi^# (5 ' 13) 

i.e. /;(A;) goes grows like k l for small k and goes to a constant for large A;. Most visibility functions one might consider 
are peaked around a typical distance which we denote by r*. The transition region between these two asymptotic 
forms of eq (5.13) occurs at At* ~ I. In general Ii(k) is a monotonically increasing function of k. 
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Deep Large- Angle Shear 



Let us define the spectral index n by 

n=k-^P 5 (k). (5.14) 

Comparing the asymptotic form of eq (5.11) with eq (5.13) we see that the Cf is well defined if n < 1 for large k and 
n > 1 — 21 for small k. It seems likely in our own universe that n < — 1 at large k while n ~ 1 for small k so there is 
no convergence problem. 

In most cases the dominant contribution to Cf is predominantly from wavenumbers k ~ l/r*. The only 
exception to this rule that one is likely to encounter in practice is when one considers small I and deep samples. This 
arises because we expect to that on large scales the power spectrum will go to an n = 1 Harrison-Zel'dovich, i.e. 

P{k) w Ak kr uz < 1. (5.15) 

The evidence for this comes only from measurements of large-angle CMBR anisotropies (Smoot et al. 1992). Studies 
of galaxy clustering have not yet reached large enough scales to see the spectrum approach n « 1. With weak lensing, 
if r,/Z > r HZ , we expect to probe larger scales than r HZ . In this case the integral of eq (5.11) will have it's a largest 
contribution from the region where n « 1 and h{k) is approximately constant, i.e. the large k limit of Ii(k) but the 
small k limit of Ps(k), and the integral will be approximated by 



C? « 367r d 




Vs. 



r* > lr HZ . (5.16) 



Since the contribution to this part of the integral is only logarithmically dominant we never expect this to be an 
extremely accurate approximation. On the other hand the this expression is not very sensitive to r HZ or the way the 
power spectrum turns over from n = 1 at small scales. The contribution to C® comes from the wavenumbers in a 
broad range from k ~ i/r* to k <~ l/r HZ rather than from a narrow range of wavenumbers near k ~ Z/r*. In our 
universe it is likely that the lo-Z shear averaged over a depth of a gigaparsecs will receive significant contribution from 
scales of a few hundred megaparsecs to gigaparsecs. 

Volume Limited Samples With Power Law Spectra 

As an example we compute the angular power spectrum of the shear in a volume limited sample of galaxies in 
the case where the density power spectrum is a power law in wavenumber, i.e. 

P(k) = Ak n . (5.17) 

For a pure power law both C? and C? have an ultra-violet (large k) divergence if n > 1 and an infrared (small k) 
divergence if n + 21 < 1. Combining the visibility function of (4.42) with the asymptotic approximation of eqs (A6) or 
(B9) one finds the approximate result 

MA-,. .("-»)(;-»?■('-■■> ( m)"" r^-v aj= 

1-21 <n<l. (5.18) 



Jl (z + i) 1 -" ~ ~ r(8-n) " U(¥W W 



324 7r 3 (11 -3n )( 6-n ) ^r (^)^ 



We illustrate that the error in this approximate result is small in fig 2. We see that for n < 1 eq (5.18) is an excellent 
approximation for all / while for smaller n there arc only significant corrections from small I. In our universe the 
density power spectrum is not an exact power law but may be approximately so at small and large scales. It is 
probable than on scales <C 100/i _1 Mpc that n w — 1 while at large scales nwl, We see that eq (5.18) diverges as 
n — > 1. This logarithmic divergence is regulated by a cutoff at small scales as explained above. 

One can argue that the asymptotic approximation of eqs (A6) and (B9) is a very good approximation when 
applied to our own universe. This follows from the fact that for deep galaxy catalogs it is the large scales which are 
relevant for the small /. At these large scales n « 1, which is where eq (5.18) is extremely accurate for all /. It is true 
that this asymptotic approximation may be off by a <~ 30% at I = 2 for extremely shallow surveys, but as we shall see 
below one needs a fairly deep survey to measure the shear. 
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Small-Angle Limit 

The small-angle shear has been studied extensively elsewhere and provides a useful check on our results so far. 
The large-^ limit of in eq (A7) is appropriate for small angle correlation functions. Taking this limit we find that 
(5.9) becomes 

C?(r,0«18**/(^yr.. , Sfi-^-^W) (5-16) 



max(i,^) 



Thus eq (4.37) becomes 

(7 2 > - <- 2 > 

and using 

j 

a "V" VI + z 

we finds that the rms shear at the redshift z is 



H ° r 2(1- ) (5.18) 



(5.19) 



3 

(1 \ ^ r ^ z^ 00 



which agrees with eq (83) of Blandford et al. (1991) and eq (2.3.6) of Kaiser (1992). Note that the polarization, p, 
used in these papers is twice the shear, 7, used here. Furthermore Blandford et al. (1991) uses units where Ho/a = 1 
while Kaiser (1992) uses units where \Ho/a = 1. Finally note that P(k) in Kaiser (1992) is defined to be a factor 
(27r) 3 larger than that used here, while Blandford et al. (1991) use the same Fourier convention as here. 

Distance-Limited Sample at Small Angles 

Now let us consider the average shear in a spatially uniform sample going to a maximum distance r max , i.e. use 
the visibility function of eq (4.42). In the small angle limit we obtain 

'HoV r°° dk {_ I 3 3/ : 

Uf ~ 107T t 1 — 

V a 

and thus 



a? „ is* l^) J^(i + - ^-J «w <».») 

'"max 

:{ 1 - :; /. 1 : ' ~ 



1 - 



TT l-OO 

— / dkkPs(k) 
a Jo 



35 V VI + 2m, 

Comparing (5.21) to (5.17) we see that the average shear in a distance limited sample is roughly 0.64. . . times the 
shear at the edge of the sample. 

Power Law Power Spectra for a Distance-Limited Sample 

Now let us also specify the power spectrum, choosing a power-law overdensity correlation function 

£(r) = (<5(x)<5(x')) - (^) 3+ " r=|x-x'| (5.22) 
which is equivalent to the power-law power spectrum 

ftw - fc " - 3< " <0 - (5 ' 23) 



Substituting this into (5.20) we obtain 



4 / , s 3+n 



/-Y f M 324,^(11 ~ 3n) (6 - n) 1 ^ ^(S - n) ^ ^ -*<»<<>. 
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which is equivalent to 2nd expression in (5.18). The reason for the divergence at n = is because the correlation 
length, r , is ill-defined for n > 0. Apart from this n-dependcnt divergence of the correlation length the /-dependence 
is correct for n < 1. For power law spectra the shear two-point correlation is given by 



C-(0, <p, <p') » C-(tf ) (cos 2{<p - <p') + (2 ^ )(4 + n) cos 2{<p + <p') 



; V 2 ; V 7 V 7 -2<n<0. (5.25) 



The correlation functions diverge for n < — 2 due to an infrared divergence. The integrals which give these correlation 
functions are formally convergent at large I for any n, however for n > —1 the contribution to the integral is dominated 
by the behaviour of the integrand at I >• l/i?. Thus for large n the results obtained are sensitive to our assumption 
of exact power law behaviour over a broad range of scales. 

One must be careful about evaluating these expressions for integer values of n. While it may appear that there 
is a singularity at n = 0, note that is zero at n = 0, and C— is smooth at n = 0. Note that these power-law spectra 
cause the correlation function be be singular at zero separation which leads to certain pathologies. For example 
is negative for all $ if n > while the coefficient of cos 2(ip + <p') in C- has a positive divergence as -d — > 0. If we 
were to regulate this divergence at small scales we would find that C« would be positive and finite at d = while the 
coefficient of cos2(<y9 + ip') in C— would go to zero. 

While the power spectrum of mass fluctuations may indeed be a power law down to very small scales, we can 
only measure shear averaged over a finite region of the sky and this averaging decreases our sensitivity to the the 
divergent density fluctuations at large scales. On small scales, the convolution of the shear pattern with a Gaussian 
beam corresponds to multiplying the angular power spectrum by exp(— l 2 a 2 ) where a is the "Gaussian width" of the 
beam, which is related to the FWHM by a f w h m = 2\/21n2<7. For a Gaussian beam we obtain 

ai q \ i a , r(2±a)r(-§)r(l-n) / flow \ 4 / r \ 3+n F ,2 + n i tf 2 1 
CW(t?)«8lV^(ll-3n)(6-n) 2 2 lFl (—-;l;- ——(5.26) 



r(3±H)r(8 - n) \ a J V2W/ V 2 ' ' 4<t 2; c7 2 +™ 

where iFi is a generalized hypergeometric function. Here we have used the notation Cx{"&) to indicate that we are 
taking the correlation of the smoothed quantities. Simpler expressions can be found for the variance of zero-lag (i.e. 
D = 0) smoothed quantities : 

r(2±2)r(-2)r(i-n) /% mra \ 4 / r \ 3+ " 1 



CW(0) w 81^(11 -3n) (6-n)/ 1+ ™ - ! ; - 
%(0, V) ¥?') w CW(0) cos 2(<^ - ^') 



r(2±2)r(8-n) V a 7 V2w7 T 2+n -2<n<0. 



(5.27) 

6. Accuracy of Shear Measurements 

When the shear is weak, i.e. when 7 <C 1, the way one goes about measuring shear is by looking for alignments 
in the apparent orientation of background galaxies. The implicit assumption is that there is no significant intrinsic 
alignment of galaxy cllipticities and any apparent alignment must be due to a coherent shear in the image. As discussed 
above this is not liable to be exactly true but is liable to be a excellent approximation. One can measure the apparent 
shape of galaxy, g, using the ellipticity tensor (Kaiser 1992) which is 

ab QaQlc-QlQl 2° ab 2 6 \sm2x 9 ~ cos2 X 9 J ( ' 

where Q 9 , Q 9 , Q 9 ab give the 0th, 1st, and 2nd moments of the galaxy brightness distribution on the sky. Here x 9 1S 
position angle of the galaxy and e 9 e [0, 1] it's the ellipticity. If the true galaxy position angles are randomly oriented 
then the expected ellipticity is 

« b ) K = {l-)f)la b (6.2) 
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in the weak lensing approximation. Here (• • -) g indicates an ensemble average over galaxy orientations, and e 2 gives 
the mean square ellipticity. Given the range of possible ellipticities the proportionality constant between ellipticity 
and shear varies by only a factor of 2 between perfectly circular galaxies, e 9 = and highly elongated galaxies, e 9 = 1. 
The variation from the expected value is 

(( e ab - ( e ab)g)( e ?d ~ ( e cd)e))e = ^e 2 S gg/ (g ac g bd + e ad e bc ) (6.3) 

where again g ac is the metric on the sphere of the sky and e a b the Levi-Civita symbol. For perfectly circular galaxies, 
e 9 = 0, the apparent position angle and shape gives the shear precisely, while for non-circular galaxies there is 
considerable variance. Typical intrinsic galaxy shapes have e 2 ~ (0.3) 2 . Note that if one makes random measurement 
errors in the shapes of galaxies then one should add this in quadrature to the deviation from the expected value in 
eq (6.3). Given the fairly large intrinsic error due to galaxy non-circularness one would have to make rather large 
measurements errors to significantly increase ones uncertainty in the determination of the shear. Roughly speaking 
one needs N ~ e 2 /[8(l — ^e 2 ) 2 7 2 ] objects to obtain a determination of the mean shear in a patch of the sky with 
S/N^ 1 if one has perfect measurements. Since typically one is looking for 7 ~ 0.01 one needs significantly more than 
100 perfectly measured galaxies or even more when one includes measurement errors. Of course the main problem in 
practice will be non-random measurement errors, i.e. errors in galaxy shapes which are correlated between different 
galaxies. 

Below we will consider estimating the shear from a uniform all-sky survey of galaxies. We will allow for an 
arbitrary weighting of different galaxy types, given by a parameter w g for each galaxy. This weighting might depend 
on the apparent shape (not orientation!) of the galaxies, the colors, or the apparent magnitude. Given this weighting 
there are a variety of different definitions of the number of galaxies in one survey, in particular one may define it with 
different powers ofw g , i.e. 

/ \ 

N gn = X /g = ( £ < Si2) ( 6 - 4 ) 



where <5-functions arc normalized such that 



J d 2 nS (2 \h-hg) = 1. (6.5) 

and (• • -) g indicates an ensemble average over galaxy positions and orientations, but not over realizations of the shear 
distribution. This is only possible because we will be assuming that the galaxy distribution and shear distribution are 
independent. We do not expect this to be exactly correct, however it is probably a fairly good approximation since the 
galaxies which sheared are so far from the galaxies associated with the mass-distribution which is doing the shearing, 
that the correlations are fairly weak. 

An estimator of the radially weighted shear is given by 

7ah(n) = 7 — r • (6.6) 

AT gl (l-Ie 2 ) 

One should, of course, not take the (5-functions as meaning the shear in concentrated in small patches of the sky, rather 
one should average this estimator over an angular patch large enough to contain many galaxies, and thereby obtain 
an estimator of the angle-averaged shear. Using the weak lensing formula, (6.2), we find that 



\E g WgS^(n-h g )j oo 

(9ab(n)) g = (lab(r g , n)) g — * = J drV(r) 7ab (r, n) - 7 afc (n) (6.7) 



where we have defined the visibility function, V(r): 
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Here we have assumed isotropy and homogeneity on the sky, so that the radial and angular distribution of galaxies 
factorizc. This shows that 7 ab (n) is an unbiased estimator of 7 oh (n). Of course one needs to know the distribution of 
galaxies in ones sample to determine V(r) and thus to know what 7 a b(n) is an unbiased estimator for. In the weak 
lensing limit we may use eq (6.3) to determine the mean square error in this estimator, which is 



6ab( r , n) - J ab (r, n)) fi cd (/> n ') ~ 7 c d( r '> «'))) = 7^2 ~, " 71 {9ac9bd + e ac e bd ) <5 (2) (n - n'). (6.9) 

(i-H 



Comparing this correlation function with C— of §4 we find the correlation function of the error in the estimator of 
eq (6.6): 

C«>(0,<p,<f/) = —&- _ cos 2(<p-<f/) lim<f(l-e-cos0) (6.10) 

4iV g i (l-I e 2j 

where the superscript "sn" refers to the shot noise from the finite number of galaxies. We may also compute the power 
spectrum for the shot noise by substituting C™ into (4.31): 

' ' "4(/ + 2)(Z-l)JV gl ^ 1 _p^- (6 - U) 

We see that the shot noise contributes equally to the scalar and pseudo-scalar part of the error. The shot noise is 
white noise, i.e. on small scales (I » 1) both C® sn and C® sn become /-independent. Substituting (6.11) into (4.27) 
we may compute the mean square error in the shear for all modes with I < L is 

<^ 2 > - ^ , e2 _, 2 X> + 1) = ^ : (L + 3)(L - 1) L > 2. (6.12) 

At small angles the mean square shot-noise in a Gaussian beam is 

SiVgl ^ _ l g 2 j JO iVgl ^ _ l g 2 j ^fwhm 

where we have used cr fwhm = 2\/21n2(7. 

Estimators for C,® and C ; ® 

A sum of unbiased estimators of some quantities is itself an unbiased estimator of the sum. Thus combining 
(6.6) with (4.12) we can construct unbiased estimators of the spherical harmonic amplitudes which describe the shear 
field 

*<•■»> -- 2 (lT2j! JVgl (i_p)t U ' fl ^'-' 

^•- a £-»)« 1 • (6 ' 14) 



+ 2)!iVgi ^_^E^^, :ah (n 9 )< h 



The superscript g emphasizes that these are estimators derived from a galaxy sample. One can construct estimators 
of the Cf and Cf using eq (4.13) 

^ S = -cr + ^TTiy E kl^T ^ s = -^ sn + wry E |<C)T- ( 6 - 15 > 

These are unbiased because the shot noise has been subtracted and differ from those in eq (4.13) since those were 
estimators derived from perfect knowledge of the shear pattern on the sky, while these include additional uncertainties 
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due to finite galaxy sampling. Of course real surveys will probably never cover the entire sky and other estimators 
must constructed to deal with finite sky coverage. 

The mean square deviation of Cf and Cf from Cf and Cf will depend on the 4-point correlation function 
of the shear as well as the clustering properties of the galaxies. We will not go into such detail here, instead we assume 
Gaussian statistics, i.e. set the reduced 4-point function to zero, and assume Poisson sampling of the shear field. These 
approximations are liable to be good ones for deep samples on large angular scales. In any case, better estimates of 
the uncertainties should be computed. With these assumptions we find 





-cf) 


(c® s - cf) 




-Cf) 


(cf £ -cf) 


w 


-Cf) 


(cf £ -cf) 



*8'2zTi(c? a + cr 9 ) 



(6.16) 



so that the 1-a fractional uncertainty is 
AC,® 



cf sn \ acv 



\ 21 + 1 V Cf 2 J Cf \ 21 + 1 y Cj 



' (6.17) 



r.2 



The 1st term in the radical represents cosmic variance which is more dominant for small I, while the second is 
uncertainty due to shot noise which is more dominant at large I. With the above estimate of the error one can define 
a total signal-to-noise: 

S V_Y- ( Cf V _^2l + l Cf 



n® j ^ Vac,® ) ^ 2 £>© 2 + £f©sn 2 

J (6.18) 

s<S y C f y _^2i + i cf 



v / ^ \ACf J ^ 2 cf 2 +Cf sn " 



which should be greater than unity in the range of I of interest if one is capable of detecting a significant shear. Note 
that these signal-to-noise ratios include cosmic variance. 

7. A Realistic Example 

Now let us estimate the angular power spectrum of the shear we might expect to find in our own universe. 
Clearly there are large uncertainties in the distribution in mass in our universe so we cannot predict precisely what 
one will see. In fact the utility of weak lensing is to measure the mass distribution. Here we give a rather simple 
empirically based model, which is meant to be optimistic about the size of the shear which one might find. We are 
also not very realistic about the data one is liable to collect. For example we assume an all-sky volume limited sample 
of galaxies, while one is really liable to obtain a partial-sky magnitude-limited sample. More detailed studies clearly 
should be done. 

On small scales the correlation function of optically selected galaxies appears to be approximately a power law 
(Peebles 1980) as was assumed for the mass in eq (5.22-27). Taking parameters. r = 5 h~ 1 Mpc and n = —1.2, which 
are close to the values obtained observationally (Humit et al. 1996) we find 

(1 \ 2 - 2 / 1 o \ °- 8 / -I o \ 0.8 

1 ~ 7f+c) (*) "•«"«">•— "{7 

C-(?9, ip, ip') w C-(tf ) (cos 2(ip - ip') + 0.5833 cos 2(ip + if')) 

/ 1 \ 2-2 / , o \ 0.8 / i o \ 0.8 

C=(0) ^ (0.0235) 2 1 - = « (0.0110) 2 z max 2 2 . (7.1) 

C=(0, <^')~C=(0) <XX2{fp - if') 
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where expressions proportional to z max 2 2 are appropriate for z max <C 1. This power law behaviour certainly does 
not persist to very large scales and it is generally believed that the spectrum turns over to something close to a 
Harrison-Zel'dovich spectrum (n = 1) at large scales. As a model of the power spectrum of density inhomogeneities 
on both small and large scales we will use 

Ps(k) - Ak A r H ^rm0.6) 

where we again use r$ = 5 /i _1 Mpc and n = —1.2 and supplement it with r HZ = 38.8/i _1 Mpc which matches the COBE 
normalization for adiabatic initial conditions. 

In fig 4 is plotted C® from this model for various survey depths, z max , and compare this with the result we would 
obtain if we used the asymptotic approximation of eq (5.9) combined with pure power law spectrum of eqs (5.22-24) 
on all scales. This figure illustrates that the approximations used to derive eq (7.1) are accurate on all angular scales 
for shallow surveys, and on small angular scales power law model for deeper surveys. In fig 5 we plot the rms shear 
contributed by angular scales larger than a given scale, as a function of the minimal scale, again for various values of 
depth. On the largest angular scales one may only find a shear of <~ 10~ 4 . We will argue below that there will be 
no significant signal for z max < 0.1 so we see that for interesting surveys one is looking for a shear > 3 x 10~ 4 . For 
such a level of shear one needs > 10 5 galaxies with good shape information to detect it. One might need many more 
if the galaxy shapes have large statistical or systematic errors. Of course as one goes deeper one needs fewer galaxies. 
One critical issue is whether one is able to make absolute measurements of the shear. If one can then one is sensitive 
to all angular scale larger than the area one surveys. Absolute measurements are possible in principle as one can use 
stellar images to correct for instrumcntally or atmospherically induced shear (Kaiser, Squires, and Broadhurst 1995). 
If for some reason one is limited to differential measurements of shear then for small survey areas one loses much of 
the signal. Except for small I the 7 rms w K rms so that fig 5 also gives the image amplification. Broadhurst, Taylor, 
and Peacock (1995) have proposed a methodology to statistically measure amplification due to gravitational lensing 
and applied it in Broadhurst (1995). 

Some will argue that we the simple mass model used here is probably an overestimate of the magnitude of 
the shear. The author considers it an optimistic estimate, but not wildly so. For non-adiabatic perturbations the 
large scale density inhomogeneities are less than assumed here. Furthermore it has long been argued that the mass 
inhomogeneities on the 8 /i _1 Mpc scale are smaller than the galaxy inhomogeneities (e.g. Henry & Arnaud 1991, White, 
Efstathiou, and Frenk 1993), so that we may also be overestimating the lensing on small scales as well. Furthermore 
the model power spectrum used here clearly has some contribution from non-linear clustering, which is not quite 
consistent with the linear evolution of density perturbations we have assumed in §5. Corrections for the non-linear 
evolution of the power spectrum will be important for deep samples, z > 0.5, on small angular scales. If galaxies are 
biased tracers of the mass distribution or if the universe is not flat then, with fixed Sp/p, the lensing shear at small 
rcdshifts scales proportionate to flo/b, where b is the bias factor. The reader may wish to apply this correction using 
his favorite values, which for many probably means reducing the predicted shear by at most a factor of 2 on small 
angular scales. 

Here is considered a volume limited survey, up to redshift z max , and assume no evolution in the number density, 
which is given by $*, and use uniform weighting of all galaxies within the survey volume, i.e. w g — 1. Using eq (6.8) 
we find 



2c ( x \ 



1 3 



= 2.6 x 10 9 ( 1 - - 1 ) (7.3) 



Zm< 

where we have taken $* = 0.012 /i 3 /Mpc 3 which is the value of this parameter in the Schechter which describes 
observed galaxies (Loveday et al. 1992). Taking e 2 sa (0.3) 2 we find 



c es„ = c ® m = 3 10 -n _([+__) ( i _ _L= ) (7.4) 



-3 



so that 



( 7 ^ L 2 ) = 6.9 x 10- 6 y/(L + 3)(L - 1) (l - -^===) '. (7.5) 
The minimum shot noise occurs when one includes only the quadrupole and at maximum depth, i.e. L = 2 and 



z = oo, which yields y (7™ 2 ) — 4.8 x 10 6 . This is extremely small, much smaller than the expected signal on these 
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scales. While going to very large redshifts may never be feasible, going to z = 1, which seems feasible, brings this up 
to 3.0 x 10~ 5 . Measurements to date have detected shear at the level of a few percent. We see that by doing very 
large area and very deep surveys one may in principle measure the shear down to a level up to 1000 times smaller. 
One may run into uncorrectable systematic errors long before one reaches this level, however it is certainly worth the 
effort get as close to the shot noise limit as possible. 

Comparison of Shot Noise With Signal 

It is interesting to compare the shot-noise "signal" with that which one might expect to obtain from density 
fluctuations. Of course density fluctuations produce no pseudo-scalar shear, so we need only compare the scalar 
component of the shear. Comparing eq (5.20) with eq (6.11) we find 



i»l 2 S e 2 V a / Jl/r k \ l ( kr n ]i 



(7.6) 



max } 



where we have used eq (A7) for the small angle limit and eq (5.20) for the volume-limited survey. This ratio is plotted 
in fig 6 and we see that for our model this signal-to- noise exceeds unity in each angular wavenumber, I, for all-sky 
surveys with z max > 0.2. Of course, as one decreases the sky coverage this signal-to-noise goes down. 

If fig 7 is plotted the fractional uncertainty in a measurement of Cf according to eq (6.18). This includes cosmic 
variance, i.e. the additional uncertainty in one's determination of the global average of quantities due to the sample 
variance associated from looking at the universe from one vantage point. From this fractional uncertainty one can 
compute the total signal-to- noise defined in eq (6.18). This signal-to-noise includes from cosmic variance, unlike the 
ratio of eq (7.6). It is perhaps most interesting to consider the ^ where it is close to unity. This leads us to consider 
shallow all-sky surveys where the contribution to from each individual / is much less than unity. From fig 7 it is 

clear that |£- should grow greater than unity for 0.1 < z m ax < 0.2. For this range of z max it is clear from fig 4 that 
the power law approximations of eqs (7.1) are a good approximation for nearly all I. We see from fig 7 that cosmic 
variance gives a negligible contribution to AC® for all I with z max this small. Furthermore it is clear that one can 
approximate the sum of eq (6.18) with an integral and take the large-Z limit of both C® and Cf sn . Combining all of 
these approximations we find that for an all-sky survey 

ce\ 2 /- 00 c® 2 /7 N 10.4 
dl I 



N®y J 2 cf 



(^) ' WJS0.2. (7.7) 



The /-integral converges relatively slowly, so that half the contribution comes from I < 10 and half from I > 10. The 
depth at which this signal-to-noise exceeds unity is not strongly dependent on the assumed amplitude of inhomo- 
geneities, for example if we were decrease the density inhomogeneities by a factor of 2 this would change value of z max 
at which ^ equals unity from 0.11 to 0.15. For fractional sky coverage, / s ky, we expect 05 fj, , although this 
depends in detail on the geometry of the region surveyed. If one makes absolute shear measurements on a small patch 
of the sky then one has limited ability to determine on which angular scale the shear is generated. The extremely 
strong dependence on z max exhibited by eq (7.7) breaks down for z max > 0.2 as one starts to sample the turnover in 
the density power spectrum. 

8. Summary and Future Development 

In this paper was presented a formalism of tensor spherical harmonics on the sphere which may be used to 
describe a weak lensing shear pattern over the entire celestial sphere. Many useful formula are presented in this 
paper but little in the form of derivations. A more complete e xposition can be found on the world wide web at 
|http: / / www-astro-thcory.fnal.gov/Pcrsonal/stebbins/WcakLcns~ . It was shown that in the all-sky shear pattern may 



be decomposed into geometrically distinct components: scalar shear and pseudo-scalar shear. In the weak lensing limit 
only scalar shear is produced by density inhomogeneities, and this is liable to be the only significant component of the 
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shear field. Pseudo-scalar shear may be produced by vector and tensor perturbations, although the amplitude is liable 
to be negligible. On small scales one may also get significant pseudo-scalar shear when the lensing becomes strong due 
to more than one mass concentration along the line-of-sight. Given that one expects negligible pseudo-scalar shear on 
the sky, it's primary use is liable to be as a gauge of measurement errors. 

In §5 the formalism was applied specifically to density inhomogeneities in an Einstein-deSitter cosmology, and 
in §6 to the study the shot noise one will obtain when trying to measure the shear with a finite number of galaxies. In 
§7 we apply the formalism derived to a model of the density perturbations in our universe. We illustrate that there is, 
in principle, a significant signal for wide area galaxy surveys even if one does not have a very deep galaxy sample. For 
example an all-sky survey of galaxies with redshifts less than 0.2 should exhibit a very significant signal of shear from 
gravitational lensing if one could reduce systematic errors below the 5 x 1CP 4 level in shear. A concerted observational 
effort should be made to see just how small a shear can be reliably measured with present technology from the ground. 
It should be noted that shallow surveys do benefit from the fact that the galaxies subtend a significantly larger angle 
than in deep surveys and thus there is less of a problem with seeing a finite pixelization. Given the large signal-to-noise 
which is present it is clear that some day such shallow large area shear maps will be made, but one may have to wait 
for large space-based optical surveys of the sky to be made. 

The estimates of the expected shear presented in this paper are rather crude in that they assumed a volume 
limited survey and did not properly take into account the evolution of nonlinear clustering. Furthermore the error 
estimates assumed Gaussian statistics. Shallow surveys even over large areas sample the nonlinear clustering regime 
which is non-Gaussian. Proper estimates of the accuracy with which one can measure the shear should include the 
proper sample variance for the non-Gaussian density field. Again further study is needed. 

Both shear and the linear polarization of light are described by a traceless symmetric rank-2 tensor on the 
celestial sphere and most of the mathematic formulae presented here can equally well be applied to a description of 
the linear polarization light on the celestial sphere. Kamionkowski, Kosowsky, and Stebbins (1996) have applied the 
tensor harmonic decomposition to a description of the cosmic microwave background radiation (CMBR). Just as with 
shear, density inhomogeneities (i.e. scalar modes) cannot, in linear theory, produce pseudo-scalar linear polarization 
of light. However in contrast to shear, vector and tensor modes can produce polarization of amplitude comparable to 
that produced by density inhomogeneities. Thus one may interpret any pseudo-scalar component of linear polarization 
of the CMBR as a direct detection of vector and/or tensor modes. This is in contrast to CMBR anisotropy where one 
cannot distinguish anisotropy induced by scalar modes from that produced by vector or tensor. 

Acknowledgements: The author would like to thank Josh Fricman, Donn MacMinn, Uros Seljak, and Nick Kaiser 
for useful conversations. Special thanks to Marc Kamionkowski who has looked over some of the formalism. This 
work was supported by the NASA grant NAG 5-2788. 

Appendix: Calculating h{x) 
In this paper we have introduced the function 
rl dy 



V¥ t fl c ,,/_J + 2 3 iI _ x\ I p ,itUl3 3 x 2 x (AJ) 



f d 

h{x) = — (1 - y)ji(xy) 
Jo V 

2<+ir(§ + i) x [1 lF2( 2 ; ~2~' 2 +l; ~ (7TT) lF2( ^~ ; 2 +l; ~T } 

where \F 2 is a generalized Hypergeometric function which can be written in terms of sines, cosines, and Sine integrals. 
Specific cases are 

-2x + x cosx + sinx + x 2 Six 



h(x) 
h{x) 



2x 2 

1 x cosx — sin x + x 2 Six 

3 2x 3 / A2 x 

-16x 3 - 3x(10 - x 2 ) cosx + 3(10 + x 2 ) sin x + 3x 4 Six 
= 24^4 

2 3x(14 + x 2 ) cosx - (42 - llx 2 ) sinx + 3x 4 Six 
"15 8~^ 
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where the sine integral function, Si, is 

f x dx 

S\(x) = I — sinx. (A3) 
Jo x 

These expressions become increasingly more complicated for large I and can involve enormous cancelation between 
terms. While symbolic algebra program do provide a generalized Hypergeometric functions, they are generally noisy 
for large arguments. Luckily one can derive accurate methods of calculating h(x) based on the asymptotic expansions 
of this function. Using the Taylor series expansion of the spherical Bessel function in the definition of we find 

(-1)" x 2n+l 

h ^ = 2 ™( 2 C + n) + 1)!! n\ (2n + T) (2n + I + 1) ' (A4) 
For large x it we can expand I\ (x) in powers of 1 jx yielding 

Between the series of (A4) for x < I and the series (A5) for x > I one can compute h{x) accurately for all values of x 
as long as / is not too large. While these series involves significant cancelation for x <~ I one can use symbolic algebra 
programs to exactly compute the series up to a given order for rational values of x and interpolate in between. 

Asymptotic Approximation and Small-Angle Limit 

For large I one finds that h(x) is well approximated by just the asymptotic form of (A5) where this is positive, 
and where (A5) is negative h(x) is negligibly small, i.e. 

where HQ is the Lorentz-Heaviside function which is unity for positive argument and zero otherwise. Taking the large 
I limit of the factorial ratios this becomes 

Il{x) *\[Il il'l) n{x ~ l) />>L (A7) 

The the fractional error of (A7) goes to zero for large I and x > I while the absolute error goes to zero for large I and 
x < I. 

Volume Average: Ii 

In this paper we have considered the mean shear from uniform distance limited in an Einstein-deSitter cosmology 
where the gravitational potential is assumed to be constant. In this case one use Ik defined in eq (5.12) with the 
visibility function of cq (4.42) to compute the mean shear. The Taylor series of the resultant function is 



^ (— 1)" x max 2 ™ + ' 

/l(-w) = 3 g 2™ (2(1 + n) + 1)!! n! (2n +"o(2n + I + 1) (2n + I + 3) Xmax = (A§) 



Note that by making /; a function of fcr max rather than just k we are using a slightly different notation than in the 
body of the paper. One may defined an asymptotic approximation just as in eq (A6) which is 



/z(2w) w (ai - -— i- + - 2 l —s U x max - - (A9) 



O '92 3/1 max 

where ai and bi are defined just as in eq (A6). The large I limit of this function is 

ij(»max) W 7 fl-~ — + n( Xmax -l) (A10) 



Z y 2Z Vy 2 x may r 2 x max ^ r 
Numerically the function // as a function of x max is not very different than the // as a function of x. 
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FIGURE CAPTIONS 



Figure 1: In the top left panel is shown a scalar pattern of gravitational deflections. Scalar deflections are parallel 
(or anti-parallel) to the direction of variation of the deflection. In the top right panel is shown a pseudo-scalar pattern 
of gravitational deflections. Pseudo-scalar deflections are ±90° from the direction of variation of the deflection. For an 
arbitrary pattern of gravitational deflections one will transform the scalar components of the deflections into pseudo- 
scalar components and vice versa by rotating the deflection direction by 90°. Plotted in the bottom left panel is a scalar 
pattern of cllipticities generated by applying the deflections in the top left panel to a square array of circles. Scalar 
shear produces ellipticity patterns which are parallel and perpendicular to the (horizontal) direction of variation of the 
cllipticities. Scalar shear produces images which arc magnified or demagnified. In the bottom right panel is shown the 
pattern of ellipticities induced by a the deflections in the top left panel. Pseudo-scalar cllipticities oscillate between 
being +45° and —45° from the direction of variation. In the weak lensing limit pseudo-scalar shear does not amplify 
the images. One will transform the scalar components of the shear pattern into pseudo-scalar components of shear 
and vice versa by rotating the object cllipticities by 45°. Such a transformation will not change the magnification or 
the displacement of the objects on the sky, so rotated scalar shear will still exhibit (de)magnification while rotated 
pseudo-scalar shear will not. While shear patterns cannot always be treated as the result of a deflection it is usually 
a good approximation to do so. 

Figure 2: Here is illustrated the meaning of the variables in the 2-point correlation function C 7 , C 7K , and C 7W . The 
thick solid line represents the geodesic connecting the two points on the celestial sphere in question, n and n', whose 
length is #. The thin solid line at n gives the orientation of the component of the shear that is being correlated, 
which is rotated ip from the direction of the geodesic. For the shear-shear correlation function, C 7 , we also need the 
orientation of the component of the shear at n' which is given by ip'. The handedness of the coordinate system is 
important C lul since it changes sign if — > ~ip. The rotation of <p> is to the right in a right-handed coordinate system. 
Note that the usual spherical polar coordinates, (9, <p) are usually left-handed on the sky since we view the celestial 
sphere from the inside. 

Figure 3: Plotted as a function of I, for a power law spectrum of density perturbations, is the ratio of the true Cf 
to that computed with the asymptotic approximation of eq (5.9). The power law index, n, is 0.9 for the crosses, for 
the circles, -1 for the triangles, and -2 for the squares. The errors continue to become smaller for larger The ratio for 
1 = 1, which contributes nothing to the shear, is not shown for n < — 1 since the true value diverges for these spectral 
indices. The asymptotic approximation is accurate except for very small /. 



Figure 4: Plotted as a function of /, for a volume limited survey in an Einstein-deSitter universe is C l obtained 

for an unbiased COBE-normalized model of adiabatic density fluctuations given by eq (7.3). The yellow, green blue, 

cyan, magenta, and red points are for a limiting survey redshift of z max =0.1, 0.2, 0.5, I, 2, and 5 respectively. The 



corresponding solid lines give the C; for the same limiting redshift but using the asymptotic approximation and 
assuming the the small scale power law index n = —1.2 on all scales. The deviation of the solid lines from the points is 
primarily due to the turnover in the spectrum at large scales, which greater for deeper surveys and smaller / because 
these probe larger comoving scales. The large I behaviour is described by eq (7.2), i.e. Cf oc l~ 1 ' 2 z max 2 - 2 , while for 
deep enough surveys, z max > 1, eq (5.16) describes the behaviour, i.e. Ci oc / with only a weak depth dependence. 

Figure 5: Plotted versus L is the rms shear contributed by modes with angular wavenumber I < L. As in fig 3 these 
are for the power spectrum model of eq (7.3) and the yellow, green blue, cyan, magenta, and red points are for a 
limiting survey redshift of z max =0.1, 0.2, 0.5, 1, 2, and 5 respectively. The the solid lines are to guide the eye. For 
I 3> 1 this also gives the rms expansion, k, which is gives the source amplification. 



Figure 6: Plotted versus I is the ratio of the shear power C l to that induced by finite galaxy sampling in a volume 
limited survey with with rms ellipticity 0.3 and galaxy density 0.012(/i/Mpc) 3 . As in figs 3&4 these are for the power 
spectrum model of eq (7.3) and the yellow, green blue, cyan, magenta, and red points are for a limiting survey redshift 
of z max =0.1, 0.2, 0.5, 1, 2, and 5 respectively. 

Figure 7: Plotted versus I is an estimate of the fractional uncertainty (1-cr) in the measurement of Cf from a volume 
limited sample with rms ellipticity 0.3 and galaxy density 0.012(/i/Mpc) 3 . This uncertainty, taken from eq (6.17), 
includes the effects of finite galaxy sampling and cosmic variance. As in figs 3 — 5 these are for the power spectrum 
model of eq (7.?). The yellow, green blue, and cyan points are for a limiting survey redshift of z max =0.1, 0.2, 0.5, and 
1. The smallest fractional uncertainty, for z max = 1 is dominated by cosmic variance for the Vs plotted. One cannot 
significantly reduce this by going to a greater depth. Less deep surveys are significantly degraded by finite galaxy 



26 



numbers. A sample with z max = 0.1 survey would yield a marginally significant result by combining all ^-modcs. One 
would not obtain a significant signal for z max significantly less than 0.1. 
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